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ABSTRACT 

We calculate evolution, collapse, explosion, and nucleosynthesis of Popula- 
tion III very-massive stars with 5OOM and 1OOOM . Presupernova evolution 
is calculated in spherical symmetry Collapse and explosion are calculated by 
a two-dimensional code, based on the bipolar jet models. We compare the re- 
sults of nucleosynthesis with the abundance patterns of intracluster matter, hot 
gases in M82, and extremely metal-poor stars in the Galactic halo. It was found 
that both 500M Q and 1OOOM models enter the region of pair-instability but 
continue to undergo core collapse. In the presupernova stage, silicon burning re- 
gions occupy a large fraction, more than 20% of the total mass. For moderately 
aspherical explosions, the patterns of nucleosynthesis match the observational 
data of both intracluster medium and M82. Our results suggest that explosions 
of Population III core-collapse very-massive stars contribute significantly to the 
chemical evolution of gases in clusters of galaxies. For Galactic halo stars, our 
[O/Fe] ratios are smaller than the observational abundances. However, our pro- 
posed scenario is naturally consistent with this outcome. The final black hole 
masses are ~ 23OM and ~ 5OOM for the 5OOM and 1OOOM models, re- 
spectively. This result may support the view that Population III very massive 
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stars are responsible for the origin of intermediate mass black holes which were 
recently reported to be discovered. 

Subject headings: nuclear reactions, nucleosynthesis, abundances - stars: evo- 
lution - stars: supernovae: general - stars: abundances - galaxies: starburst - 
galaxies: intergalactic medium 

1. Introduction 

One of the most interesting challenges in astronomy is to investigate the mass and 
properties of first generation "Population III (Pop III)" stars, and how various elements 
have been synthesized in the early universe. Just after the Big Bang these elements were 
mostly only H, He and a small amount of light elements (Li, Be, B, etc). Heavier elements, 
such as C, O, Ne, Mg, Si and Fe, were synthesized during the evolution of later generation 
stars, and massive stars exploded as supernovae (SNe), releasing heavy elements into space. 

Stars end their lives differently depending on their initial masses M . Here the Pop III 
stars are assumed to undergo too little mass loss to affect the later core evolution. Then 
the fates of Pop III stars are summarized as follows. Those stars lighter than 8M form 
white dwarfs. Those with 8M - 130M Q undergo ONe-Fe core collapse at a last stage 
of their evolution leaving neutron stars or black holes. Some of these stars explode as 
the core-collapse supernovae. Stars with 13OM - 3OOM undergo electron-positron pair 
creation instability during oxygen burning, releasing more energy by nuclear burning than 
the gravitational binding energy of the whole star, and hence these stars disrupt completely 
as the pair-instability supernovae (PISN). Stars with 3OOM - 1O 5 M also enter into the 
pair-instability region but continue to collapse. Fryer, Woosley, & Heger (2001) calculated 
evolution of 260 and 300 M stars and obtained the result that a 260 M star ends up as a 
PISN and a 300 M star collapsed. Stars over ~ 1O 5 M collapse owing to general relativistic 
instability before reaching the main-sequence. The core collapse SNe (Type II, lb and Ic 
SNe) release mainly a-elements such as O, Mg, Si and Ca and some Fe-peak elements as 
well. 

It has been suggested that the initial mass function (IMF) of Pop III first stars may 
be different from the present one - that more massive stars existed in the early universe 
(e.g., Nakamura & Umemura 1999; Abel, Bryan, & Norman 2000; Bromm, Coppi, & Larson 
2002; Omukai & Palla 2003). Some authors (e.g., Wasserburg and Qian 2000; Qian, Sargent, 
Wasserburg 2002; Qian and Wasserburg 2002; & Yoshida et al. 2004) argued that existence 
of very massive stars (VMSs) in the early universe is consistent with abundance data of Ly« 
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systems. Numerical simulations by, e.g., Bromm & Loeb (2004), indicate that the maximum 
mass of Pop III stars to be formed will be ~ 300M Q - 5OOM . Omukai & Palla (2003), 
however, point out that under certain conditions VMSs much heavier than 300M Q can be 
formed in the zero-metallicity environment. Tan & McKee (2004) calculated star formation 
by taking rotation and disk structure and concluded that first stars should be much more 
massive than 3OM . Another scenario for the formation of VMSs for any metallicity has 
been presented by Ebisuzaki et al. (2001) and Portegies Zwart et al. (1999, 2004a, 2004b), 
where VMSs are formed by merging of less massive stars in the environment of very dense 
star clusters. 

In the present paper, we call the stars with M ^ 1O 5 M "Super- Massive Stars (SMSs)", 
and the stars with M = 13OM - 1O 5 M "Very-Massive Stars (VMSs)". Among "VMSs" 
we define M > 3OOM stars as "Core- Collapse Very-Massive Stars (CVMSs)", in order to 
clarify the distinction between the PISN mass range and the core-collapse range. Here we 
focus on CVMSs, and deal with 500 and 1000 M models. 

Such CVMSs might have released a large amount of heavy elements into space by mass 
loss and/or supernova explosions, and they might have significantly contributed to the early 
galactic chemical evolution, if they were the source of reionization of intergalactic H and 
He (e.g., Gnedin & Ostriker 1997; Venkatesan, Tumlinson, & Shull 2003). The reionization 
of intergalactic He has traditionally been attributed to quasars. However, according to the 
results of the the Wilkinson Microwave Anisotropy Probe (WMAP) observation in 2003, 
reionization in the universe took place as early as 0.2-0.3 billion years after the Big Bang 
(redshift z ^ 20)(Kogut et al. 2003). Then these Pop III CVMSs might provide a better 
alternative channel which could operate at redshifts higher than what is assumed for quasars 
(Bromm, Kudritzki, & Loeb 2001). 

The question of whether CVMSs (~ 3OOM — 1O 5 M ) actually existed is of great impor- 
tance, for instance, to understand the origin of intermediate mass black holes (IMBHs) (~ 
5 x 1O( 2-4 )M ). Stellar mass black holes (~ 1OM ) are formed as the central compact 
remnants of ordinary massive (25 - 13OM ) stars at the end of their evolution, while super- 
massive black holes(SMBHs) (~ 10 5 - 1O 9 M ) are now known to exist in the center of almost 
all galaxies (e.g., Kormendy & Richstone 1995; Bender 2004). IMBHs have not been found 
until recently. However, there is a strong possibility that some IMBHs have been, indeed, 
found (e.g., Barth, Green, & Ho 2005 for most recent review). Matsumoto et al. (2001) 
reported possible identification of a ^ 7OOM black hole in M82, by using Chandra data. 
As to formation of SMBHs there are several scenarios (e.g., Rees 2002, 2003). SMBHs may 
be formed directly from supermassive halos of dark matter (e.g., Marchant & Shapiro 1980; 
Bromm and Loeb 2003). Ebisuzaki et al. (2001) suggested a scenario where IMBHs grow to 
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a SMBH by merging and swallowing of many of these objects. If CVMSs actually existed, 
they could be considered as natural progenitors of IMBHs. 

Motivated by these backgrounds, here we calculate the evolution, collapse, explosion, 
and nucleosynthesis of Pop III CVMSs (over 300M Q ). These stars are expected to form black 
holes directly at the end of evolution. It has not been known yet if they will explode as SNe. 
However, if the star is rotating the whole star will not become a black hole at once, but it 
is expected to form an accretion disk around the central remnant (e.g. Shibata & Shapiro 
2002). After forming an accretion disk, jet-like explosions may occur by extracting energy 
from the accretion disk and/or the black hole itself (Fryer et al. 2001; MacFadyen et al. 
2001; Maeda & Nomoto 2003). Therefore, in our current explosion and collapse calculations 
we adopt a two-dimensional approach including accretion along the equatorial direction and 
jets toward the polar direction. 

We compare our results of nucleosynthesis with the observed abundance data of intra- 
cluster medium (ICM), intergalactic medium (IGM), gases in the central part of M82, and 
extremely metal-poor (EMP) stars in the Galactic halo. Since it is very difficult to observe 
directly the explosions of the first generation stars due to the large distance (redshift z ^ 20), 
currently comparison of the kind carried out in this study will offer a powerful method to 
support the existence of such very massive stars. 

After describing the basic methods adopted and the assumptions made for our models 
in section 2, the results are presented and discussed in section 3, and they are compared 
with observations in section 4. Further discussion and concluding remarks are given in the 
last section, section 5. 

2. Methods, Assumptions & Models 

We calculate evolution, core collapse, explosion, and nucleosynthesis of very-massive 
stars with 5OOM and 1000M o . As mentioned in Introduction, such massive stars may be 
formed in a metal-free environment. We start our evolutionary calculations by assuming 
that the stars have 500 and 1OOOM , with zero-metallicity, on the pre-main sequence. As 
our starting approximation we neglect radiative mass loss due to zero metallicity (Kudritzki 
2000). Ibrahim et al. (1981), Baraffe et al. (2001), and Nomoto et al. (2003) showed that 
pulsational mass loss is not so effective for metal-free stars, so we also neglect the pulsa- 
tional mass loss. To calculate presupernova evolution we adopt the stellar evolution code 
constructed by Umeda & Nomoto (2002) based on the Henyey method. This code is devel- 
oped from the codes constructed by Nomoto & Hashimoto (1988), and Umeda et al. (1999). 
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The nuclear reaction network for calculating nucleosynthesis and energy generation at each 
stage of the evolution is developed by Hix & Thielemann (1996). We include 51 isotopes up 
to Si until He burning ends, and 240 up to Ge afterwards. Our evolutionary calculations 
are carried out from the pre-main sequence up to the iron-core collapse where the central 
density reaches as high as 2 x 10 10 gcm~ 3 . When the temperature reaches 5 x 10 9 K, where 
"nuclear statistical equilibrium" (NSE hereafter) is realized, the abundance of each isotope 
is determined for a given set of density, temperature, and Y c . Here Y e is the number of 
electrons per nucleon, defined as: 

n = £§*. (i) 

where Zi is the atomic number, Ai is the mass number, and Xi is the mass fraction of species 
i. Y c , as well as density and temperature, is a key quantity to determine the abundance of 
each element. We assume NSE at logT (K) > 9.7. 

The explosion mechanism of core-collapse supernovae is not well understood. Moreover, 
we do not know beforehand how very massive stars over 300M Q explode due to strong 
gravitation even when they are rotating. In this study, therefore, instead of going into 
the problem of whether such massive stars actually explode, we investigate the conditions 
required for these stars to explode, by exploring several situations with various models. 
For the explosion in the hydrodynamical simulation, we adopt the two-dimensional (2D) 
Newtonian hydrodynamical code constructed by Maeda & Nomoto (2003) and Maeda (2004). 
This code adopts the Eulerian coordinate and solves Euler equations based on Roe's scheme 
(Hachisu et al. 1992, 1994). Previously 2D simulations of jet-induced supernova explosions 
have been carried out by many authors for ordinary massive stars with ~ 25M Q to 40M Q 
(Khokhlov et al. 1999; MacFadyen et al. 2001; Nagataki 2000; Maeda et al. 2002; Maeda & 
Nomoto 2003; Maeda et al. 2006). However, there have been no such detailed calculations 
for CVMSs with M ^ 5OOM . Because temperatures are so high at the explosion, the 
pressure is radiation-dominated, and hence we use the equation of state for the radiation 
and electron-positron pairs with the adiabatic index of 4/3. 

The explosion models to be explored are summarized in Tables 1 and 2. Jets, which are 
supposed to be injected from the accretion disk, are considered to be the energy source of the 
explosion. We first choose the initial black hole mass, and the outer matter accretes toward 
the central object. Because our hydrodynamical code includes gravitational force, the final 
black hole mass and the ejected mass are determined as the results of the calculations for a 
set of given parameters (Maeda & Nomoto 2003). One of the purposes of this study is to 
explore the condition for the stars to explode when jets are injected. As typical cases we 
choose the initial black hole mass, M BH o, at 100M o for the 1000M Q star and 5OM for the 
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5OOM star. However, in order to investigate the dependence of results on this parameter, 
we also explore larger or smaller values of M B ho- 

We adopt the spherical polar coordinate with the number of meshes set to 150 x 60 for 
five models (higher resolution; A-l, A-2, B-l, B-2, and B-3), and 100 x 30 for the rest (lower 
resoluton). The latter models are chosen so that we can search for parameter dependence 
quickly. We describe our results mainly based on those of higher resolution models. However, 
we include lower resolution models to explore the detailed dependence on parameters. We 
also calculate Model B-2 with lower resolution, in order to numerically compare the two 
resolutions. The difference of the quantities we give in Tables 3, 4, 6, and 7 between the 
two resolutons is around 10 %. 

At the beginning of hydrodynamical simulations, the region Mbho < M r is mapped onto 
the computational domain. The central part (M r < Mbho) is displaced by a point mass. 
The inner boundary of the simulations is set at the radius Ro (see Tables 1 and 2). In the 
computational domain, we assume that the effect of rotation is negligible. This assumption 
applies if the specific angular momentum j 17 = j/10 17 cm 2 s" 1 in the progenitor star is in 
the range 6.3 < jn « 45, where the lower and upper limits correspond to the conditions 
that the disk forms beyond the schwarzshild radius and well below the inner boundary of 
our computational domain. If jn ~ 6.3, which is favorable in order to make an efficiently 
accreting disk (Narayan, Piran, & Kumer 2001), then the rotational force is at most a few 
percent of the gravitational force at the inner boundary. 

For the jet injection, we choose various values for the parameter, 9 jet, the angle from 
the polar axis. The jet is injected into the direction of < 9 < 9j Ct . At the direction 
9 > 9 j et , the inner boundary is treated as follows: it is set to be transmitted (absorbed; i.e., 
vanishing radial gradient of all variables) or reflected if the material just above the boundary 
has negative (accreting) or positive sign. These boundary conditions are used in previous 
studies on jet-induced supernova explosions (Khokhlov et al. 1999; Macfadyen, Woosley, 
& Heger 2001, Maeda & Nomoto 2003). It should be noted that by using the transmitted 
boundary condition for the accretion case and neglecting pressure and rotational support 
below the boundary, we may overestimate the accretion rate. 

The energy and mass injected by the jets per unit time are connected with the properties 
of accreting matter as (Maeda & Nomoto 2003): 

Ejet = eM^ecC 2 = thermal + (^Pjct^ 2 cl>jet-4jct (2) 
Mjet = I^M acc = PjctVjctAjet (3) 

where M acc is the accretion rate, Ej et the injected energy per unit time, Mj et the mass spouted 
per unit time, e the energy transformation efficiency, /i the mass fraction of the jets to the 
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accreted matter, pj ct the jet density, i>j e t the jet velocity, and A$ et the area over which the 
jet is spouted, respectively. We treat e and ji as free parameters to be varied to explore the 
explosion energy. 

We consider two cases for the form of the injected energy. One is that almost all the 
energy of the jet is given as kinetic energy (Case A). The other is that almost all the energy 
is given as thermal energy (Case B). We introduce a parameter -F therma i defined as: 



i.e., the ratio of thermal energy in the jet to the total jet energy per unit time. This parameter 
is set to 0.01 for Case A, and to 0.9 or 0.95 for Case B. By using equations 2, 3 and 4 we 
obtain the jet velocity: 



Equations 2,3, and 4 give a complete set of jet properties at the inner boundary. In 
Case A the jet carries most energy toward the polar (jet-injected) direction, and hence the 
model is highly non-spherical. For this case, we set the ratio between the two parameters, 
e//i, to 0.1. Because we perform a Newtonian calculation, the larger we set e the larger we 
need to set p, so that the jet material does not exceed or approach the speed of light. For 
Case B, on the other hand, the models, due to the dominant thermal energy, become more 
spherical, because thermal motion is random and non-directional. For this case, we set larger 
e values for the same /i compared with Case A. This means that the jet is something like a 
hot bubble. 

One of our primary purposes is to investigate how much heavy elements are synthesized 
and ejected by the explosion. Therefore we stop the calculations when all of the following 
conditions are satisfied: 

1 M a cc decreases enough below O.1M s _1 , i.e., ^ O.O2M s -1 . 

2 Total explosion energy E to t becomes much larger than the absolute gravitational bind- 
ing energy |£ grav |, E tot Z 10|£ grav |. 

These criteria mean that accretion has almost stopped. 

3 The maximum temperature of the matter decreases below 8 x 10 8 K. 

This means that explosive nucleosynthesis no longer occurs at such low temperatures. 



'thermal 



(4) 



thermal 





(5) 



Under these criteria, calculations sometimes end before jets reach the stellar surface. 
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Maeda & Nomoto (2003) used helium stars (the hydrogen envelope is removed by mass 
loss) as the initial models, and carried out the calculations for about 100 s until the jet reaches 
the stellar surface and the expansion becomes homologous. In contrast the radii of the stars 
we use here are in the order of lO 2 i? (~ 10 13 cm) because they have the hydrogen-rich 
envelope. Because we investigate the first generation stars, the mass loss will not be effective 
due to the metal-free environment. Therefore, it is reasonable to examine a hydrogen-rich 
star rather than a He star. 

We calculated explosive nucleosynthesis by using temporal histories of density and tem- 
perature stored during hydrodynamical calculations. The reaction network we use includes 
280 isotopes up to 79 Br. At high temperatures, T 9 = T/10 9 K > 5, NSE is realized. We use 
the "NSE" code (Hix & Thielemann 1996) for T 9 > 6. 

The jet matter should be included in the ejected matter and we need to calculate its 
nucleosynthesis. We do not know which of the accreted matter is injected as jets, and so 
the final chemical composition is uncertain. However, Pruet et al. (2004) carried out nucle- 
osynthesis of disk wind for various Y c values. MacFadyen & Woosley (1999) and MacFadyen 
(2003) also considered disk wind. Based on these works in this study we make the following 
assumptions: 

1 Because the jet matter is injected through the inner region (from the accretion disk), 
it should have experienced high temperatures at which NSE is realized (T 9 > 5). 

2 The jet matter expands adiabatically after it is injected (i.e., entropy is conserved). 

3 The accreted matter is mostly accreted while the accretion rate M acc is of the order 
of 10 - 1O 2 M s _1 . It is likely that p varies depending on when the jet material is 
injected. Therefore s (entropy density, oc ^-) is likely to vary as well. 

4 For the value of Y c we assume 0.48 < Y c < 0.52. 

Based on these assumptions we start the calculation of nucleosynthesis of the jet matter 
from Tg = 6, using the historical temperature and density data of the first test particle of the 
jet (injected at the first stage of the explosion) in Model A-2 (history A), and the changed 
entropy data (double (history B) or triple (history C) the density at the same temperature). 
In other words, we use three p—T histories (history A:(p(t), T(t)), history B:(2p(t), T(t)), and 
history C:(3p(t), T(t))), where the set (p(t),T(t) ) is given by the hydrodynamic simulations. 
Y e is parameterized at 0.48, 0.49, 0.50, 0.51, and 0.52. We calculate 15 patterns and average 
these results to the first approximation. Y c and entropy of the jet material can change when 
it is ejected. Therefore here we consider combination of jets with different values of these 
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parameters. These assumptions still include large uncertainty, but our aim is just to roughly 
estimate the amount of 56 Ni. The larger the mass of the jet is, the larger we expect the 
uncertainties to be. 



3. Results 
3.1. Presupernova Evolution 

3.1.1. Evolutionary Tracks 

Figure 1 shows the evolutionary tracks of the central density - temperature relation for 
5OOM and 1OOOM stars. We also plot, for comparison, the track of a 300M Q star, which 
results in the pair-instability supernova. Generally, more massive stars have higher entropies 
(lower densities) at the same temperatures (i.e., at the same burning stage). Although each 
star passes through the region of electron-positron pair-instability, both 5OOM and 1OOOM 
stars proceed to iron-core collapse (Fe-decomposition region in Figure 1), unlike the 3OOM 
star. The 5OOM and 1OOOM stars do not become pair-instability supernovae though they 
pass through the pair-instability region, because the energy released at this stage is less than 
the gravitational binding energy of the star (Rakavy et al. 1967; Bond et al. 1984; Glatzel 
et al. 1985; Woosley 1986). 

3.1.2. Presupernova Model 

Figure 2 shows the presupernova chemical composition for the 5OOM and 1OOOM 
stars. In the region labeled as "NSE region", NSE is realized. For this region, we calculate 
the evolutional changes in terms of (Y e , p, T) to obtain the NSE abundances. One can see the 
onion-like structure from the center to the surface - the iron-core, silicon layer, oxygen layer, 
helium layer, and hydrogen-rich layer. Here we define the iron-core as the region where the 
mass fraction of Si is less than 10%. The iron-core occupies up to 13OM of mass from the 
center for the 5OOM star, and 25OM for the 1OOOM star. For both cases, they occupy a 
quarter of the total mass. This fraction is much larger than that in ordinary massive stars. 
For example, in a 25M star, the iron-core is about 1.6M (Umeda et al. 1999), less than 
10% of the total mass. The reason is the difference of the density and temperature structure. 

Figure 3 shows the density and temperature structure of the two stars just before the 
explosion (when the central density reaches 10 10 g cm" 3 ), which is compared with the 25M 
model. The density and temperature gradients for the 5OOM and 1OOOM stars are smaller 
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than those of the 25M© star, and hence the regions with high temperature and high den- 
sity are larger. Then the fraction of the iron-core is larger. The large drop of density at 
M r /M total ~ 0.5 for the 5OOM and 1000M Q stars in Figure 3 corresponds to the boundary 
between the oxygen and helium layer. 

3.2. Explosion Hydrodynamics 

We describe the results of explosion hydrodynamics in this subsection and nucleosyn- 
thesis in the next subsection, showing several figures. All figures are based on the results of 
high resolution models, except for Figure 4. 

3.2.1. Explosion Energy and Ejected Mass 

In Tables 3 and 4 we summarize for each model the total explosion energy, final black 
hole mass, and mass of the jets. The total explosion energy is of the order of 10 54 erg for 
most cases, except that in Model A-4 it is of the order of 10 53 erg. Model A-4 is almost at the 
border between the 'successful' and 'failed' explosions. Actually, we also try to calculate the 
case which has the same parameters as Model A-4 except that e and fi are half the values of 
A-4 (see Model F-l in Table 5), but in this model the jet promptly falls back to the central 
remnant after it is injected, and hence the explosion fails. In this model, the total energy is 
still negative and the stellar matter moves toward the central remnant more than 200 s after 
the beginning of the accretion. The absolute value of gravitational binding energy over the 
region outside of the central 100M o core is as high as 10 55 ergs for the 1OOOM model. In 
this case energy injection is too weak for the jet to proceed outward. 

Table 5 summarizes the models in which explosion ends up as 'failure'. Figure 4 shows 
the models in which the explosion either occurs or not, depending on the two parameters 9j et 
and e. One can see that the minimum e needed for the successful explosion becomes higher 
if 0jet is larger, as in Model A-2 and Case B (most of the injected energy is given as thermal 
energy). Actually, explosion energies tend to be lower in such models than those in Case 
A models with 9- ]et = 15°. In Model A-2 and Case B models, the injected energy tends to 
diffuse into the direction apart from the polar direction, and so the jet is weak even for the 
polar direction. 

Figure 5 shows how the jet is propagating through the star by plotting the density 
structures. One can see that the jet is strongly propagating in the polar direction for A-l, 
while for B-l the jet is broadened toward the side directions due to the random heat motion. 
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The final black hole mass and ejected mass are also important. For the 500 M Q and the 
1OOOM models, these values are ~ 23OM and 5OOM , respectively. 

3.2.2. Direction- Dependent Features 

Figure 6 shows the regions where the matter will be accreted onto the central black hole. 
In these figures we can see the extent of asphericity and the amount of accreting matter. 
These panels show the initial positions (just before the explosion) of the accreted matter. 
In Model A-2 in which #j et is twice (30°) the other models (15°) and Case B, asphericity is 
weakened to some extent and the amount of the accreting matter is less compared with the 
models with 9- ]Ct = 15° and the same energy transformation efficiency e. For models with 
#jct = 15°, the stellar matter toward the direction 9 > 9j et almost accretes up to 500M Q 
(even a part of the helium layer) on the mass coordinate. On the other hand, in Model A-2 
and Case B models a large amount of matter inner than the 5OOM core is ejected. At the 
same time, if asphericity becomes weaker the threshold efficiency for the successful explosion 
becomes more strict- that is, larger e is needed. This is because more energy diffuses toward 
the equatorial direction. 

In Figure 7 we show the maximum temperatures which each mesh reaches and the 
densities at the maximum temperatures for the z (polar)-direction, 9 = 15°, 9 = 45°. We 
may pay special attention to the maximum T 9 for the 9 = 15° and 9 = 45° direction. For some 
models (e.g., A-l in the left panel) with 9- )et = 45°, the maximum T 9 does not appear within 
the range of these graphs, because the matter which can experience such high temperatures 
is in the inner region, and hence such matter all accretes in these models. However, for the 
other cases, the inner matter is ejected and the explosive nucleosynthesis occurs even for the 
9 = 45° direction. 



3.3. Explosive Nucleosynthesis 

When a shock arrives the shocked region is compressed and heated, drastically raising 
the density and temperature, and then the explosive nucleosynthesis occurs. The products of 
this event are characterized by the peak temperature. We first summarize the main products 
at different peak temperatures and then describe the results of the calculations. 
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3.3.1. Explosive Burning and Products 

If the peak temperature T pea k exceeds 5 x 10 9 K, NSE is realized. In such regions 
'complete silicon burning' occurs and then Fe-group elements (such as Mn, Co, Fe, Ni) are 
produced. The main product is 56 Ni, which eventually decays into 56 Fe. 

In the complete silicon burning region, at lower density for a given temperature the 
reaction rate decreases, and the number of free-particles may exceed the NSE value. Or, if 
the initial temperature is higher, free-particles become more abundant because in NSE the 
number of these particles is a high-powered function of temperature. This situation is called 
'a-rich freezeout', and it tends to produce the Fe-group elements and nuclei to the high-Z 
side of the peak (e.g. Thielemann, Nomoto, & Hashimoto 1996; Arnett 1996). 

If 4 x 10 9 K < Tp Ca k < 5 x 10 9 K, incomplete Si burning occurs. In such regions, Si is not 
all converted into the Fe-group elements but remains or is converted to the elements such as 
32 S, 36 Ar, and 40 Ca. 

If 3 x 10 9 K < Tp Ca k < 4 x 10 9 K, explosive oxygen burning occurs, which produces 28 Si 
and 32 S, while the original 16 composition stays the same. 

If 2 x 10 9 K < Tp Ca k < 3 x 10 9 K, explosive carbon burning occurs, which produces 20 Ne 
and 24 Mg. The original 12 C remains because the burning does not proceed during such a 
short time scale. 

If Tp e ak < 2 x 10 9 K, almost no explosive burning occurs, and so the original chemical 
composition realized during the hydrostatical burning phase is conserved. 

3. 3. 2. Direction- Dependent Features 

As typical interesting cases, Figures 8 - 10 show the distribution of elements after the 
explosive nucleosynthesis for Models A-l, A-2 and B-l, respectively. In each figure the upper 
left panel shows the Fe-group elements in the polar direction. The upper right panel shows 
the a - elements in the polar direction, the lower left panel shows the a - elements at 9 = 15°, 
and the lower right panel shows the a - elements at 9 = 45°. In each model, complete silicon 
burning region shows strong a-rich freezeout. The upper panel (polar direction) in each 
figure shows that 56 Ni is synthesized dominantly up to 4OOM from the center. Compared 
with Figure 2 which shows the chemical composition just before the explosion, one can see 
that oxygen is consumed in the region with 350 - 400 M Q . 

For directions 9 = 15° and 9 = 45° the silicon and oxygen layers considerably accrete, 
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and even a part of the helium layer accretes for 9 = 45° for Model A-l. On the other hand, 
in Models A-2 and B-l the complete silicon burning region still remains for 9 = 15° and the 
oxygen layer remains for 9 = 45°. This is because the shock is diffused to the equatorial 
directions more than Model A-l. 

Figure 11 shows the mass fractions of 56 Ni and a-elements for each 9 integrated over 
the radial direction. These figures clearly show how much matter is ejected - for example, 
the ejected mass of the oxygen layer by seeing the mass fraction of 16 0. For the 9- ]et = 15° 
models in Case A (e.g., A-l) there is no ejected matter except helium and hydrogen in the 
directions 9 > 45° (see also Figure 8). However, for Case B and the 9j et = 30° models of Case 
A (e.g., A-2 or B-l) the 56 Ni synthesized region and a-element rich region are broadened to 
around 9 = 30° and 80°. 



3.3.3. Composition of Jet Material 

Figure 12 shows [X/Fe] (top panel) and mass fractions (bottom panels) for the Fe-peak 
elements as a function of Y e . Note that we assume the temperature of the jet material reaches 
higher than 5 x 10 9 K, and therefore it consists mostly of the Fe-group elements and 4 He. 

Peculiar features are seen particularly when Y e < 0.5. Co, Cu, Ni and Zn are dra- 
matically abundant relative to Fe (500-1000 times larger than the solar values) as shown in 
Figure 12. When Y e < 0.5 the mass fraction of synthesized 56 Ni is very small, less than 10% 
for Y c = 0.49 and less than 0.1% for Y c = 0.48. Then a large amount of neutron-rich nuclei, 
such as 58 Ni, 60 Ni, and 59 Ni (which decays into 59 Co), 63 Zn (which decays into 63 Cu), and 
64 Zn, are synthesized. In these situations neutron-rich 64 Zn is directly synthesized rather 
than 64 Ge, which decays into 64 Zn. The rise of [Cr/Fe] and [Mn/Fe] for Y c < 0.5 is mainly 
due to the small fraction of 56 Ni rather than the increase of the fractions of Cr and Mn. 

On the other hand, for Y c > 0.5 most of the products are 56 Ni and 4 He, similar to the 
case where Y e = 0.5. The main effect of Y e larger than 0.5 is the existence of free protons. 

As our first step for the treatment of the jet material, Figure 13 shows the abundance 
pattern of jet material averaged over 15 patterns (5 Y c values x 3 entropy values). [Zn/Fe] 
and [Ni/Fe] are larger than the solar values due to the effects of small Y e regions, while 
[Cr/Fe] and [Mn/Fe] are smaller due to the effects of large Y c regions. The averaged mass 
fraction of 56 Ni is about 40%. We multiply the mass fraction of each nucleus by the jet mass 
and add it to the total abundance pattern. 
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3.4. Ionization Rates, Heavy Element Yield, and Ionization Efficiency 

The suggestion that VMSs are responsible for the reionization of HI and Hel is not a 
new one (e.g., Gnedin & Ostriker 1997). Bromm, Kudritzki, & Loeb (2001) calculated the 
stellar atmosphere models for Pop III main-sequence CVMSs of 300 - 1000 M and obtained 
the effective temperatures of log T eS (K) ~ 5.05, which are higher than log T eff (K) ~ 4.81 of 
Pop I stars with the same mass and slightly higher than log T cff (K) = 4.85 - 5.0 for Pop III 
15 - 90 M® stars (Tumlinson & Shull 2000). Thanks to the high effective temperature, Pop 
III CVMSs give high production rates of ionizing radiations ~ 1.6 x 10 48 photons s _1 M _1 
for H I ionization, 1.1 x 10 48 photons s _1 M & ~ 1 for He I ionization, and 3.8 x 10 47 photons 
s _1 M Q ~ l for He II ionization. These numbers correspond to ~ 16, 14 and 75 times higher 
than the corresponding values with a Salpeter IMF (see Bromm et al. 2001), and therefore, 
they are sufficient for completely reionizing IGM. 

Daigne et al. (2004) estimate the efficiency of supplying UV photons and chemical 
enrichment of IGM simultaneously. These authors suggest that the IMF that essentially 
formes less than 100 M is favorable. However, this conclusion is due to the assumption 
that all CVMSs collapse entirely to a black hole. Venkatesan & Truran (2003) considered 
the relation between the reionizing radiation and metal enrichment of IGM, using stellar 
atmosphere models and model yields available at that time. For the model yields, they 
assumed no metal ejection by stars of ~ 30 — 130M Q and also M > 3OOM . Following their 
argument, here we compute the reionization efficiency for our CVMSs using model yields in 
the present work. 

Adopting the mass of heavy elements ejected by our 1000 M star model, M z ~ 5OM , 
the conversion efficiency (i] Lyc ) of energy produced in the HI ionizing radiation divided by 
the energy produced in the rest mass of metals (Mzc 2 ) is r^yc ~ 0.05. (We used Eq. 1 of 
Venkatesan & Truran 2003). Here we use the timescale of t ms = 2 x 10 6 years for the 1000 M 
Pop III star. With these values, the number of ionizing photons per baryon in the universe 
generated in association with the IGM metallicity Ziqm ~ 10~ 4 , obtained for our model, is 
N-Lyc/Ny, ~ 150. (We used Eq. 2 of Venkatesan & Truran 2003). Note that this value well 
exceeds the value required for reionization of inter galactic hydrogen, 1 < N^/Ny, ^ 10 
(see Somerville et al. 2003). Therefore, contrary to the earlier results, our conclusion is that 
CVMSs can contribute significantly to reionization of IGM in the early epochs. 
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4. Integrated Abundance Patterns and Comparison with Observations 

The abundance pattern, the mass ratio of each element to be compared with obser- 
vations, is determined by integrating the distributions over the entire ejecta regions (both 
radial and 6 directions). It is the mass ratio of each ejected element. Tables 6 and 7 show the 
ejected masses of some isotopes excluding the jet materials, and Table 8 shows the masses 
of all the isotopes including the jet materials. 

4.1. Abundance Patterns without Jet Materials 

Tables 6 and 7 show the ejected mass of 56 Ni (which decays into 56 Fe), excluding the 
jet material. Masses of 16 and 28 Si are also shown as representative a-elements to see the 
abundance and ratios of these elements. In models for Case A (except Model A-4), the ratios 
of the ejected masses of these elements to their progenitor mass are rather small, compared 
with those ratios in ordinary massive stars such as a 25M star. The typical ejected 56 Ni 
mass in the 25M star is ~ 0.1 M Q (Maeda & Nomoto 2003). In the models with jet = 15°, 
asphericity is so strong that it is only toward small 9 direction where 56 Ni and Fe-group 
elements are synthesized and ejected. On the other hand, in models for Case B and Model 
A-4, these masses are much larger than the other models. The ejected 56 Ni mass is about 5 - 
1OM . If this kind of supernova occurs, it is very bright in its tail because the heating source 
of a supernova is 7-rays from radioactive decays of 56 Ni ~^ 56 Co ^ 56 Fe. However, it is very 
difficult to observe directly the explosions of first generation stars by present observational 
devices since they are very distant (z ^ 20). 

4.2. Total Abundance Pattern and Comparison with Observational Data 

4-2.1. Intracluster Matter and Hot Gas in M82 

Figure 14 and Figure 15 show the total abundance pattern for each model of higher res- 
olution and lower resolution models, respectively, which is compared with the observational 
data of intracluster medium (ICM) and M82. Note that the following discussions do not 
depend on the resolution. In these figures the abundance data for the ICM gas are shown 
with the bars (Baumgartner et al. 2005; Peterson et al. 2003), while the pentagons show 
the data for the gas of the central region of M82 (Ranalli et al. 2005; see also Tsuru et al. 
1997). These data show that (1) the ratio [O/Fe] is smaller than the solar value, (2) [Ne/Fe] 
is about the solar value, and (3) the intermediate-mass a-elements such as Mg, Si, S exhibit 
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oversolar abundance - that is, [Mg/Fe], [Si/Fe] and [S/Fe] ~ 0.5 (Origlia et al. 2004; Ranalli 
et al. 2005). 

Note that these data are not explained by standard Type II SN nucleosynthesis models. 
Also, if underabundance of [O/Fe] is due to the contribution of Type la SNe, other a- 
elements such as Si and S should be also underabundant. Loewenstein (2001) suggested the 
contribution of Pop III hypernovae (supernovae of ordinary massive stars such as 25M , with 
the explosion energy of at least ~ 10 times larger than normal supernovae; e.g., Nomoto et 
al. 2003) to the enrichment of ICM, in order to explain low [O/Fe] and high [Si/Fe], using 
the hypernovae models by Nakamura et al. (2001) and Heger et al. (2001). Umeda et al. 
(2002) also discussed this feature, but they predicted smaller [Ne/Fe] and [Mg/Fe] than the 
data given by Ranalli et al. (2005). 

Here we compare nucleosynthesis calculations of our CVMS models with these observa- 
tional data. The results are summarized as follows. For our Case B models we obtain the 
abundance pattern generally close to the observations of both ICM and M82 - for instance, 
the underabundance of [O/Fe] and [Ne/Fe], and the over solar values of [Mg/Fe], [Si/Fe] and 
[S/Fe]. On the other hand, all Case A models result in very underabundant values of [a/Fe] 
& —1 because the mass fraction of 56 Ni synthesized in the jet matter is much larger than 
that synthesized in the matter which does not accrete. What is more, contribution by the jet 
material is dominant in such models and the uncertainty is very large. The yields of PISNe 
(Heger & Woosley 2002; Umeda & Nomoto 2002) are [O/Fe] ~ [Mg/Fe], [Ne/Fe] ~ 0, and 
[Si/Fe] and [S/Fe] ~ 1.0, not consistent with these data. Therefore the yields of our Case B 
models for CVMSs can explain these abundance patterns better than those of PISNe. 



4-2.2. Intergalactic Medium 

The abundances in the intergalactic medium (IGM) at high redshift also provide im- 
portant information on the early chemical evolution of the universe. Many researchers have 
attempted to measure the metallicity of IGM at high redshift (Songaila & Cowie 1996; 
Songaila 2001; Schaye et al. 2003). Aguirre et al. (2004) observed the abundances of C 
and Si in IGM at redshift 1.5 ^ z ^ 4.5, argued that Si and C have the same origin, and 
obtained [C/Si] ~ —0.77. This value is considerably lower than the yields by Population 
III ordinary massive stars (M < 4OM ) (see also Heger & Woosley 2002; Chieffi & Limongi 
2002; Umeda & Nomoto 2002). 

Matteucci & Calura (2005) discussed whether this ratio could be reproduced with their 
chemical evolution models and obtained [C/Si] = —2.0 — 1.7 by including the contribution 
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of Pop III stars over 1OOM . This is too small to be compatible with the observed value. 
Thus they concluded that the contribution of VMSs could not be large. However, they 
adopted the yields of PISNe (13O-3OOM ) only (Umeda & Nomoto 2002; Heger & Woosley 
2002). PISNe enrich much more Si than C. 

Although there is a similar feature between PISNe and CVMSs in that [C/Fe] < 0, 
[Si/Fe] > and thus [C/Si] < 0, [C/Si] from CVMSs is not so extreme as that by PISNe. With 
our yields of 1000M Q CVMS, [C/Si] ~ —0.86 — 0.68 (including both high and low resolution 
models of 1000M Q ), which is more than 10 times larger than the results by Matteucci & 
Calura (2005). Our values are compatible with the observed value (Aguirre et al. 2004), and 
they show that the contribution of CVMS to the IGM enrichment is significant. 

4.2.3. Extremely Metal-Poor Stars 

Figure 16 compares the yields of our models with extremely metal-poor (EMP) stars in 
the Galactic halo (the data by Cayrel et al. 2004). The result is that our Case B CVMS 
models mostly agree with these Galactic halo star data for both a-elements and iron-peak 
elements. [Mg/Fe] in the EMP stars is oversolar for a wide range of metallicity. [Cr/Fe] 
and [Mn/Fe] are small while [Co/Fe] and [Zn/Fe] are large. Cr (produced as 48 Ti) and Mn 
(produced as 55 Co) are mainly produced in the incomplete silicon burning region, while Co 
(produced as 59 Cu) and Zn (produced as 64 Ge) are mainly produced in the complete silicon 
burning region. We note that the aspherical models for ordinary massive stars with 25M and 
4OM in Maeda & Nomoto (2003) are also consistent with EMP star's abundance patterns. 
Umeda & Nomoto (2003) obtained similar results with spherical models by introducing a 
mixing and fallback scenario. 

It has been reported that [O/Fe] is generally oversolar for EMP stars, which does not 
agree with our models. However, there are little data for [O/Fe] at [Fe/H] ^ —3 and the 
uncertainties involved in the non-local thermal equilibrium (NLTE) effects and 3D effects 
may be too large to make conclusive statements. Therefore, to answer the question of 
whether metal-free CVMSs could contribute to the enrichment at [Fe/H] < —3, we will need 
more accurate observational data of [O/Fe]. 
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5. Summary & Discussion 
5.1. Summary 

We first calculated the evolution of Pop III CVMSs with M = 5OOM and 1OOOM 
from the pre-main-sequence through the collapse with spherical symmetry. These CVMSs 
are thought not to explode if they undergo spherical collapse (Fryer et al. 2001). We assumed 
that these stars explode in a form of bipolar jets, and explored the required constraints. The 
results of our nucleosynthesis calculations were used to examine their contribution to the 
chemical evolution of galaxies. Our major findings are: 

1 The region which experiences explosive silicon burning to produce iron-peak elements 
is more than 20% of the total mass, much larger than those of ordinary massive stars 
such as a 25M star. Note that for the metal-free 25M star model, this fraction is 
less than 10% (Umeda & Nomoto 2002). This is because for the 5OOM and 1000M Q 
models the density and temperature distributions are much flatter than those of 25M 
stars. 

2 Typical explosion energy is of the order of 10 54 erg for 1000M Q models for the parameter 
ranges in this study. 

3 Black hole masses are ~ 5OOM for the 1000M o star models. Note that such a black 
hole mass is very similar to those of IMBHs, e.g., a claimed ~ 7OOM black hole in 
M82. It is quite possible that CVMSs could be the progenitors of IMBHs. 

4 Nucleosynthesis yields of CVMS have similar patterns of [a/Fe] to the observed abun- 
dance patterns of both ICM and gases of the central region of M82 if the contribution 
of the jet is small (Case B). Specifically, for Case B small ratios of [O/Fe] and [Ne/Fe] 
combined with large [Mg/Fe], [Si/Fe] and [S/Fe] (i.e. large [(Mg, Si, S)/0]) are gen- 
erally more consistent with these observational data than those of hypernovae and 
PISNe. 

5 For IGM, [C/Si] of our CVMS models is compatible with that of IGM at high redshift 
(z — 5), which is sufficiently higher than those of PISNe. 

6 For Fe-peak elements, the main feature of the yields of our Case B CVMSs is that 
[Cr/Fe] and [Mn/Fe] are small while [Co/Fe] and [Zn/Fe] are large. This is consistent 
with the observed ratios in the extremely metal-poor (EMP) stars. The oversolar ratios 
of some a-elements, such as [Mg/Fe] and [Si/Fe], are also consistent with EMP stars. 
Our CVMS models do not agree with the oversolar [O/Fe] of EMP stars. However, 
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more data of [O/Fe] in EMP stars will be needed in order to see whether CVMSs 
can contribute to the early galactic chemical evolution. In this sense [O/Fe] would be 
important to discriminate between different models. 

5.2. Discussion 

5.2.1. Mass Accretion and Mass Loss 

It was pointed out (Omukai & Palla 2003; Tan & McKee 2004) that after a protostar 
starts shining as a main-sequence star the accretion still continues. In our current study, as 
a starting point the effect of accretion on mass growth during the presupernova evolution 
is not included. In our next more realistic models such accretion will be included in the 
evolutionary calculations also. However, Omukai & Palla (2003) find that when the protostar 
simulation of very massive stars is carried out properly with time-dependent accretion rates, 
the rates generally decrease toward the end of the protostar era and after the onset of the 
main-sequence stellar phase. 

Also, it was pointed out (Maeder & Meynet 2004) that mass loss will not be negligible 
even for zero-metallicity stars when they are rotating, and hence mass loss also will be 
included in the next step of our models. However, we expect that our major conclusions as 
summarized above are still valid, at least qualitatively. Somewhat more massive CVMSs, 
however, may be needed to obtain the same mass black holes if mass loss is significant. 

5.2.2. Reionization and Chemical Enrichment 

For our CVMSs, the timescale of evolution from the zero-age main sequence to core- 
collapse is ~ 2 x 10 6 years - only 1/3 - I/f0 as long as for ordinary massive stars (13 - 25 
Mq). So if these CVMSs were formed as the first generation stars they would be the first 
contributor to reionize and enrich the universe (Omukai & Palla 2003; Tumlinson & Shull 
2000; Bromm et al. 2001; Schaerer 2002). 

Concerning the existence of VMSs, it was proposed (e.g., see Wasserburg, & Qian 2000; 
Qian et al. 2002; Qian & Wasserburg 2002) that the prompt inventory involving VMSs 
produced the elements from C to the Fe group, in order to explain the observed jump in 
the abundances of heavy r-process elements at [Fe/H] ~ —3, and also that while VMSs 
themselves produced no heavy r-elements, these stars dominated chemical evolution earlier 
at [Fe/H] < —3. Some others (e.g., Venkatesan & Truran 2003; Tumlinson, Venkatesan, & 
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Shull 2004) argue that various observational data on reionization, the microwave background, 
the metal enrichment of the high redshift IGM, etc., indicate that the IMF of the first stars 
need not necessarily have been biased toward high masses. In what follows, we revisit this 
issue on the basis of our present models. 

In section §3.4 we estimated the ionization efficiency of our CVMSs. It was found that 
the number of ionizing photons per baryon in the universe, generated in association with the 
IGM metallicity Zigm ~ 10~ 4 , is A^L yc /A^ ~ 150, and so CVMSs can contribute significantly 
to reionization of IGM in the early epoch. Here we emphasize that our current result for 
CVMS is contributed from the mass range with ~ 300 — 1OOOM , and hence the PISN 
(pair instability supernova) range is not included. On the other hand, Venkatesan & Truran 
(2003) give N Lyc /N h ~ 10 for Z lGM ~ 10" 4 for the mass range ~ 100-1000M© which reflects 
the large contribution of PISNe to metal enrichment. Daigne et al. (2004) also considered 
reionization and chemical enrichment of IGM simultaneously and reached similar conclusion 
that VMSs are not necessary. However, in their models CVMSs do not explode. Here we 
may note that less massive Pop III stars (^ 100M Q ) can also produce the amount of ionizing 
photons per metal similar to CVMS (Venkatesan & Truran 2003; Tumlinson et al. 2004). 

The relation between reionization and metal enrichment of IGM becomes clearer if we 
solve the equation for N^yc/N^ (Eqn. 2 of Venkatesan & Truran 2003) for a given value of 
Zigm- For a 1000 M & star, Z/Z & ~ 10 -3 ' 4 and 10 -4 ' 4 for the required number of ionizing 
photons per baryon 10 and 1, respectively. This is about one order of magnitude smaller than 
the case for the mass range 100 — 1000M Q (mainly contributed by PISNe). The difference 
between CVMSs and PISNe is larger if we consider the enrichment of iron. The 260 M & 
PISN of Heger & Woosley (2002) gives Z Fe /Z Fc<Q ~ 10~ 2 - lO" 3 , while our 5OOM and 
1OOOM star gives ~ 10~ 3 - 2 - 10~ 4 - 2 . 

A main critique against the existence of PISNe comes from the fact that we do not 
see the abundance patterns of PISNe in EMP stars (Umeda & Nomoto 2002; Tumlinson et 
al. 2004). The EMP abundances are indeed suggested to be accounted for by hypernovae 
or faint supernovae of less massive stars of ^ 100M Q (Umeda & Nomoto 2003). However, 
the apparent lack of evidence of VMSs by no means contradicts the existence of CVMSs 
at earlier epochs, if the majority of first stars in the earlier epoch has masses ^ 3OOM . 
First, PISNe from stars of ^ 300M© will be just a minor fraction in such a case, explaining 
the lack of the signature of PISNe. Second, Z/Z Q expected from our CVMSs is smaller 
than PISNe. Namely, the metal enrichment by CVMSs might be finished before ordinary 
core-collapse SNe become dominant. Note that the abundance, especially of oxygen, in EMP 
stars and IGM is different. Here we have shown that the yields of our CVMSs can reproduce 
the abundance of IGM (Section 4). Therefore, it would be worthwhile studying a scenario 
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where CVMSs are first formed in pregalactic mini halos, and then subsequently ordinary 
core-collapse SNe took place in the galactic halo. 

We mentioned that the early universe would have been contaminated too much if there 
were many PISNe. This conclusion is not affected when the calculation of PISNe consid- 
ers the effect of rotation and asymmetric explosion (see Stringfellow & Woosley 1988). In 
PISNe with rotation, all matter is ejected with nothing left as in non-rotating models, and 
the total amount of matter ejected is almost independent of the geometry of the explosion, 
whether spherical or not. The explosive nucleosynthesis itself and hence the resulting exact 
composition of the ejected matter change with rotation and under the consequent asymmet- 
ric environment, but that does not change our conclusion that too many of PISNe result 
in the overabundance of heavy elements in the early universe which contradicts with the 
observation. 



5.2.3. Initial Mass Function 

Tumlinson et al. (2004) raised two problems associated with top heavy IMF of the first 
stars, (a) If stars are all M <; 300M Q , no metals are released from Pop III stars to trigger 
the transition from the first stars to present IMF star formation, and (b) no mechanism has 
been proposed for forming stars more massive than ~ 3OOM without forming PISNs. 

Here we discuss how our CVMS models could resolve these apparent problems. Con- 
cerning (a), our present CVMS models do eject metals (though less than PISNe), leading 
to metal enrichment of IGM. In this connection, note that the existing literature concerning 
the effects of VMSs (in Pop III IMF) on reionization, etc., includes only contribution by 
PISNs but not those heavier because they assumed that heavier stars do not explode, and 
hence make no contribution. However, we emphasize the importance of the explosion of 
these heavier stars (> 500M o ). 

Concerning (b), collisions and merging of ordinary massive stars in very dense clusters 
are expected to lead to the formation of more massive stars (Ebisuzaki et al. 2001; Portegies 
Zwart et al. 1999, 2004a, 2004b) which can easily lead to CVMSs (k, 3OOM ) without or 
with only minor fraction of stars responsible for PISNe (^ 300M Q ), and hence the problem 
in question can disappear. Specifically, in this scenario PISN stars will have no time to 
explode before merging into heavier stars when the timescale of the PISN star evolution is 
longer than merging timescale. Also, even in the case of single star formation (no merging) 
there is yet no reason to exclude a possibility of the first star IMF with the minimum mass 
of ~ 300M m . 
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As to the question of how CVMSs are formed, the first generation stars are generally 
thought to have formed in low-mass halos with the virial temperature T vir < 10 4 K. Then 
the upper limit to the mass of first stars may be ~ 3OOM (e.g., Bromm & Loeb 2004). 
However, Oh & Haiman (2002) investigated halos with higher mass, with T vir > 10 4 K. The 
evolution of these high mass halos, e.g., of ~ 1O 9 M , is found to be quite different from the 
low mass case, and the degree of fragmentation of the gas is still highly uncertain. Therefore, 
it appears that whether more massive stars can be directly formed is still an open question. 
However, regardless of the feasibility of direct formation, it has been emphasized by, e.g., 
Ebisuzaki et al. (2001), and Portegies Zwart et al. (1999, 2004a, 2004b), that CVMSs will 
be formed easily by merging of less massive stars in very dense star clusters, and hence there 
appears to be essentially no problem for CVMS formation. 

5.2.4. Black Hole Mass 

As to the existence of IMBHs, currently extensive effort is under way to try to detect 
them in nearby galaxies (see, e.g., Barth et al. 2005 for recent review). Already several 
of these objects have been identified, mostly in dwarf elliptical galaxies, but some in spiral 
galaxies, e.g., M33, IC 342, Pox 52, NGC 4395 and several galaxies in the Sloan Digital Sky 
Survey (SDSS) sample. For example, the black hole mass obtained for M33 is less than ~ 
3OOOM but larger than the mass of a stellar mass black hole. The black hole mass obtained 
for NGC 4395 is ~ 10 4 - 1O 5 M , but its spectra are unlike NLSI (narrow line Seyfert I) - a 
class of AGN which tends to have small mass. We expect more of these IMBH candidates, 
with better mass measurement, to be identified in the very near future. 

As an example of possibly more recently formed IMBHs, Matsumoto et al. (2001) re- 
ported the possible discovery of a ^ 7OOM black hole in M82 as an ULX (ultra-luminous 
X-ray source). Since an ULX was first detected in 1989 by the Einstein Observatory (Fab- 
biano 1989), many of these objects have been discovered. Possible scenarios for formation 
of IMBHs associated with ULXs are speculated in a recent article by Krolik (2004). Col- 
bert & Mushotzky (1999) first suggested that these luminous objects are indeed IMBHs, 
because their luminosity is super Eddington for stellar mass black holes if spherical accre- 
tion is adopted. That may not be necessary if beaming, etc., is assumed. However, ULXs 
may be heterogeneous, and at least some of these objects may well prove to be IMBHs. 
The prospect is bright because various observations in multifreuency bands can distinguish 
between different interpretations. 
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Table 1. The higher resoluton models. 



Models 


progenitor(M ) 


M BH o(M ) 


i? (km) 


9jct° 


e 


H 


-^thermal 


Ujet/C 


A-l 


1000 


100 


1.5 x 10 4 


15 


0.01 


0.1 


0.01 


0.45 


A-2 


1000 


100 


1.5 x 10 4 


30 


0.01 


0.1 


0.01 


0.45 


B-l 


1000 


100 


1.5 x 10 4 


15 


0.01 


0.005 


0.95 


0.45 


B-2 


1000 


100 


1.5 x 10 4 


15 


0.005 


0.0025 


0.95 


0.45 


B-3 


500 


50 


1.1 x 10 4 


15 


0.01 


0.005 


0.95 


0.45 



Note. - initial black hole mass M B ho; the radius at the inner boundary of the simulations 
Ro, jet injected angle #j e t, energy tranformation efficiency e, the mass fraction in which accreted 
matter is ejected as jet /x, and the jet velocity normalized by the speed of light t>j Ct . 



Table 2. The lower resolution models. 



Models 


progenitor(M ) 


Mbho(M ) 


i? (km) 




e 




-^thermal 


Ujet/C 


A-3 


1000 


100 


1.5 x 10 4 


15 


0.005 


0.05 


0.01 


0.45 


A-4 


1000 


100 


1.5 x 10 4 


15 


0.002 


0.02 


0.01 


0.45 


A-5 


1000 


50 


8.4 x 10 3 


15 


0.005 


0.05 


0.01 


0.45 


A-6 


1000 


200 


2.7 x 10 4 


15 


0.01 


0.1 


0.01 


0.45 


A-7 


500 


50 


8.4 x 10 3 


15 


0.01 


0.1 


0.01 


0.45 


B-2 


1000 


100 


1.5 x 10 4 


15 


0.005 


0.0025 


0.95 


0.45 


B-4 


1000 


100 


1.5 x 10 4 


15 


0.01 


0.02 


0.9 


0.32 


B-5 


1000 


100 


1.5 x 10 4 


15 


0.005 


0.01 


0.9 


0.32 
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Table 3. Caracteristics of explosion models in Table 1. 



Models 


E (erg) 


M BH (M ) 


M eO (M ) 


M jet (M ) 


M e (M ) 


A-l 


6.7 x 10 54 


5.0 x 10 2 


4.6 x 10 2 


44 


5.0 x 10 2 


A-2 


2.2 x 10 54 


4.4 x 10 2 


5.2 x 10 2 


38 


5.6 x 10 2 


B-1 


4.9 x 10 54 


4.6 x 10 2 


5.4 x 10 2 


1.8 


5.4 x 10 2 


B-2 


1.6 x 10 54 


4.8 x 10 2 


5.2 x 10 2 


0.94 


5.2 x 10 2 


B-3 


2.9 x 10 54 


2.3 x 10 2 


2.7 x 10 2 


0.90 


2.7 x 10 2 



Note. — Explosion energy E, final black hole mass M B h, ejected mass 
excluding jet material M e0 , mass of jet Mj et , and total ejected mass M e . 

Table 4. Same as Table 3, but for the models in Table 2. 



Models £(erg) M BH (M ) M eO (M ) M jet (M ) M e (M ) 





A-3 


4.0 x 10 54 


5.3 x 10 2 


4.5 x 10 2 


23 


4.7 x 10 2 






A-4 


5.2 x 10 53 


5.6 x 10 2 


4.3 x 10 2 


9.5 


4.4 x 10 2 






A-5 


3.0 x 10 54 


5.1 x 10 2 


4.5 x 10 2 


45 


4.9 x 10 2 






A-6 


8.2 x 10 54 


5.1 x 10 2 


4.6 x 10 2 


34 


4.9 x 10 2 






A-7 


6.1 x 10 54 


2.4 x 10 2 


2.4 x 10 2 


21 


2.6 x 10 2 






B-2 


1.9 x 10 54 


4.8 x 10 2 


5.2 x 10 2 


0.95 


5.2 x 10 2 






B-4 


6.4 x 10 54 


4.6 x 10 2 


5.3 x 10 2 


7.4 


5.4 x 10 2 






B-5 


1.8 x 10 54 


4.7 x 10 2 


5.3 x 10 2 


3.7 


5.3 x 10 2 




Table 5. The models in which explosion does not occur. 


Models 


progenitor(M ) 


M BHO (M ) 




n 


-^thermal 




F-l 




1000 


100 


15 0.001 


0.01 


0.01 


0.45 


F-2 




1000 


100 


30 0.005 


0.05 


0.01 


0.45 


F-3 




1000 


100 


45 0.01 


0.1 


0.01 


0.45 


F-4 




1000 


100 


15 0.003 


0.0015 


0.95 


0.45 
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Table 6. Ejected mass of 56 Ni, 16 0, 28 Si, excluding jet material for higher resolution 

models. 



Models 


M( 56 Ni) (M ) 


M( 16 0) (M ) 


M( 28 Si) (M ) 


A-l 


1.5 


4.5 


0.69 


A-2 


12 


18 


4.3 


B-l 


9.3 


23 


6.1 


B-2 


4.7 


24 


4.0 


B-3 


8.9 


10 


4.8 



Table 7. Same as Table 6, but for lower resolution models. 



Models 


M( 56 Ni) (M ) 


M( 16 0) (M ) 


M( 28 Si) (M ) 


A-3 


0.4 


2.1 


0.24 


A-4 


0.12 


1.2 


0.080 


A-5 


0.36 


2.4 


0.23 


A-6 


1.4 


5.1 


0.77 


A-7 


1.1 


1.8 


0.41 


B-2 


4.4 


22 


3.8 


B-4 


7.0 


24 


5.3 


B-5 


5.8 


25 


4.6 
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Table 8. Nucleosynthesis products (M ) of models A-1, B-1, B-2, and B-3. 



ivioacis 


A 1 


R 1 


R 9 
Jj-Z 


-D-O 




A 1 


R 1 

LJ- 1 


R 9 


Q 

1 > ■ ' 


n 


1.68E-13 


6.94E-12 


5.19E-13 


5.65E-14 


B0 Cr 


1.04E-03 


1.75E-04 


9.71E-05 


1.29E-04 


p 


1.90E+02 


1.94E+02 


1.94E+02 


9.59E+01 


51 Cr 


1.75E-05 


8.55E-07 


7.88E-07 


6.45E-07 


d 


2.78E-15 


3.65E-15 


2.84E-15 


4.00E-15 


52 Cr 


3.81E-05 


1.61E-06 


8.16E-07 


9.77E-07 


3 He 


2.12E-04 


2.12E-04 


2.12E-04 


2.02E-05 


53 Cr 


6.53E-09 


4.55E-10 


2.15E-10 


2.28E-10 


4 He 


2.72E+02 


2.85E+02 


2.83E+02 


1.43E+02 


54 Cr 


4.17E-10 


1.30E-10 


5.06E-11 


8.69E-11 


6 Li 


1.09E-18 


9.64E-18 


1.20E-17 


6.99E-18 


55 Cr 


7.20E-12 


1.41E-10 


4.98E-11 


5.59E-11 


7 Li 


2.77E-14 


1.36E-13 


1.52E-13 


3.47E-14 


48 Mn 


6.84E-05 


2.99E-06 


1.37E-06 


1.40E-06 


7 Be 


3.77E-09 


3.86E-09 


3.89E-09 


1.86E-09 


49 Mn 


3.54E-04 


1.12E-04 


2.88E-05 


6.88E-05 


9 Be 


6.40E-21 


3.83E-19 


3.10E-20 


8.71E-21 


50 Mn 


2.30E-04 


6.94E-05 


2.08E-05 


3.61E-05 


8 B 


2.93E-14 


3.05E-14 


3.05E-14 


4.54E-13 


51 Mn 


1.16E-03 


2.68E-04 


9.65E-05 


1.56E-04 


10 B 


7.04E-15 


8.02E-14 


9.11E-14 


2.96E-14 


52 Mn 


1.09E-04 


9.84E-06 


7.71E-06 


9.95E-06 


n B 


1.08E-15 


1.22E-14 


1.20E-14 


5.68E-15 


53 Mn 


1.66E-04 


2.00E-05 


1.37E-05 


2.87E-05 


n C 


1.74E-12 


1.05E-11 


4.28E-12 


3.14E-12 


54 Mn 


3.20E-07 


1.38E-08 


6.94E-09 


7.63E-09 


12 C 


8.11E-01 


3.56E+00 


3.68E+00 


7.61E-01 


55 Mn 


1.73E-08 


9.05E-10 


4.36E-10 


5.97E-10 


13 C 


6.07E-09 


2.23E-08 


2.30E-08 


4.58E-08 


56 Mn 


8.12E-12 


1.03E-10 


4.28E-11 


8.74E-11 


13 N 


1.16E-10 


3.20E-10 


6.14E-10 


2.31E-08 


57 Mn 


7.48E-12 


1.08E-10 


4.22E-11 


7.60E-11 


14 N 


1.16E-05 


3.49E-05 


2.77E-05 


7.15E-06 


50 Fe 


2.01E-04 


1.09E-05 


3.63E-06 


4.36E-06 


1B N 


3.13E-09 


1.30E-07 


3.93E-08 


1.40E-07 


51 Fe 


4.36E-04 


5.11E-05 


9.14E-06 


1.84E-05 


14 o 


2.54E-05 


9.14E-04 


1.86E-04 


1.41E-04 


52 Fe 


4.31E-02 


9.34E-02 


5.22E-02 


8.61E-02 


15 Q 


2.18E-06 


7.72E-06 


4.25E-06 


5.12E-06 


53 Fe 


1.36E-03 


3.34E-03 


2.00E-03 


3.61E-03 


16 


4.55E+00 


2.31E+01 


2.37E+01 


1.03E+01 


54 Fe 


2.26E-03 


5.19E-03 


4.99E-03 


1.11E-02 


17 o 


1.74E-08 


1.29E-07 


5.67E-08 


3.57E-08 


55 Fe 


3.96E-05 


3.52E-06 


2.16E-06 


1.06E-05 


18 o 


2.50E-07 


5.92E-07 


9.19E-07 


6.93E-11 


56 Fe 


5.40E-05 


2.73E-06 


1.33E-06 


6.56E-06 


17 F 


1.24E-08 


2.86E-09 


1.58E-09 


7.01E-10 


57 Fe 


6.49E-08 


2.94E-09 


1.45E-09 


1.63E-09 


18 p 


8.20E-08 


8.77E-06 


1.25E-06 


7.24E-07 


58 Fe 


6.00E-09 


4.07E-10 


1.87E-10 


2.09E-10 


19 p 


6.80E-11 


5.56E-08 


1.37E-09 


1.43E-08 


59 Fe 


6.49E-12 


1.16E-10 


5.11E-11 


6.46E-11 


18 Ne 


6.92E-06 


3.10E-07 


1.48E-07 


4.68E-07 


60 Fe 


7.62E-12 


1.31E-10 


4.37E-11 


7.54E-11 


19 Ne 


2.91E-07 


1.07E-06 


2.77E-08 


5.09E-07 


61 Fe 


3.56E-12 


1.20E-10 


2.07E-11 


7.08E-11 


20 Ne 


7.83E-01 


5.18E+00 


4.82E+00 


1.58E+00 


51 Co 


2.89E-10 


1.19E-11 


5.95E-12 


5.71E-12 


21 Ne 


1.23E-05 


2.56E-05 


4.48E-05 


2.35E-05 


52 Co 


1.22E-05 


4.82E-07 


2.27E-07 


2.35E-07 


22 Ne 


1.34E-05 


2.27E-05 


4.64E-05 


1.50E-06 


53 Co 


1.56E-04 


1.39E-05 


4.15E-06 


6.23E-06 


21 Na 


4.98E-07 


2.84E-06 


6.22E-07 


1.44E-06 


54 Co 


3.07E-04 


6.25E-05 


1.14E-05 


3.25E-05 


22 Na 


1.64E-06 


1.13E-05 


1.53E-05 


3.55E-05 


55 Co 


1.37E-03 


3.30E-03 


1.79E-03 


4.02E-03 


23 Na 


4.18E-05 


3.46E-04 


4.14E-04 


5.56E-04 


56 Co 


1.72E-04 


1.89E-05 


1.32E-05 


1.99E-05 


22 Mg 


2.37E-04 


1.58E-05 


6.22E-06 


1.26E-05 


57 Co 


5.33E-04 


2.71E-05 


1.60E-05 


1.43E-05 


23 Mg 


2.14E-04 


7.84E-04 


9.12E-04 


1.04E-03 


58 Co 


7.77E-06 


3.29E-07 


1.66E-07 


1.60E-07 


24 Mg 


7.05E-01 


4.55E+00 


4.14E+00 


1.43E+00 


59 Co 


3.23E-05 


1.37E-06 


6.89E-07 


6.58E-07 


25 Mg 


1.47E-04 


6.76E-04 


7.04E-04 


2.88E-04 


60 Co 


4.68E-09 


4.20E-10 


1.80E-10 


2.39E-10 


26 Mg 


3.65E-05 


1.70E-04 


1.50E-04 


4.11E-05 


61 Co 


1.82E-11 


2.24E-10 


8.42E-11 


1.02E-10 


27 Mg 


1.54E-11 


9.35E-10 


1.85E-10 


3.38E-11 


62 Co 


6.29E-12 


1.32E-10 


5.11E-11 


8.78E-11 


25 Al 


2.97E-05 


1.10E-04 


1.25E-05 


2.79E-04 


54 Ni 


5.26E-05 


1.42E-06 


5.54E-07 


6.72E-07 


26 Al 


3.88E-06 


1.72E-05 


1.42E-05 


2.54E-05 


55 Ni 


3.16E-04 


1.57E-05 


4.70E-06 


5.77E-06 


27 A1 


5.97E-04 


7.06E-03 


4.95E-03 


5.78E-03 


56 Ni 


2.04E+01 


1.01E+01 


5.07E+00 


9.32E+00 


28 Al 


1.98E-08 


5.58E-07 


1.88E-07 


9.17E-08 


57 Ni 


3.80E-01 


1.68E-01 


8.32E-02 


1.30E-01 


29 Al 


1.42E-12 


7.85E-11 


1.31E-11 


1.24E-11 


58 Ni 


7.66E+00 


3.83E-01 


1.93E-01 


2.03E-01 
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Table 8 — Continued 



ivioacis 


A 1 
A- 1 


R 1 


R 


-D-O 




A 1 
A-1 


R 1 




J3-Z 


R *3 
D-O 


26 Si 


5.63E-04 


7.80E-05 


2.16E-05 


5.95E-05 


59 Ni 


1.47E-01 


6.79E- 


03 


3.76E-03 


3.47E-03 


27 Si 


2.62E-04 


1.21E-03 


1.06E-03 


1.14E-03 


60 Ni 


2.31E+00 


9.74E- 


02 


5.21E-02 


4.69E-02 


28 Si 


6.92E-01 


6.19E+00 


4.04E+00 


4.81E+00 


61 Ni 


2.29E-03 


9.65E- 


05 


4.97E-05 


4.65E-05 


29 Si 


3.07E-04 


2.26E-03 


1.94E-03 


3.03E-03 


62 Ni 


2.13E-04 


9.00E- 


06 


4.64E-06 


4.33E-06 


30 Si 


1.66E-05 


2.38E-04 


7.63E-05 


5.81E-04 


63 Ni 


3.05E-08 


1.40E- 


09 


6.71E-10 


6.68E-10 


31 Si 


4.78E-11 


8.74E-09 


6.55E-10 


3.94E-09 


64 Ni 


6.81E-10 


2.02E- 


10 


9.08E-11 


1.11E-10 


32 Si 


5.46E-13 


1.14E-11 


5.42E-12 


3.31E-12 


65 Ni 


6.02E-12 


1.32E- 


10 


4.23E-11 


7.07E-11 


27p 


1.93E-06 


8.17E-08 


4.11E-08 


3.99E-08 


66 Ni 


7.20E-12 


2.09E- 


10 


5.35E-11 


8.03E-11 


28 p 


7.79E-05 


8.57E-06 


2.43E-06 


3.40E-06 


56 Cu 


1.23E-09 


5.07E- 


11 


2.53E-11 


2.45E-11 


29 p 


1.04E-04 


4.08E-04 


6.05E-05 


1.35E-03 


57 Cu 


1.29E-03 


5.67E- 


05 


2.78E-05 


2.87E-05 


30 p 


1.26E-04 


3.25E-04 


1.75E-04 


2.79E-04 


58 Cu 


3.60E-02 


4.92E- 


02 


1.41E-02 


4.46E-02 


31p 


8.43E-05 


4.22E-04 


2.66E-04 


7.31E-04 


59 Cu 


4.34E-02 


3.67E- 


02 


1.15E-02 


3.46E-02 


32 p 


1.37E-09 


5.33E-09 


3.76E-09 


1.15E-08 


60 Cu 


6.17E-03 


3.54E- 


03 


2.70E-03 


1.71E-03 


33 p 


4.73E-10 


1.16E-09 


1.25E-09 


4.37E-09 


61 Cu 


6.53E-02 


2.80E- 


03 


1.55E-03 


1.35E-03 


34 p 


1.76E-12 


2.02E-11 


1.41E-11 


2.85E-12 


62 Cu 


2.45E-03 


1.04E- 


04 


5.26E-05 


5.00E-05 


30 g 


1.05E-03 


1.58E-04 


3.32E-05 


7.83E-05 


63 Cu 


1.81E-03 


7.66E- 


05 


4.00E-05 


3.69E-05 


31 g 


2.60E-04 


2.94E-04 


9.71E-05 


1.88E-04 


64 Cu 


1.85E-06 


7.84E- 


08 


3.96E-08 


3.78E-08 


32 g 


3.55E-01 


3.17E+00 


2.00E+00 


2.68E+00 


65 Cu 


7.68E-07 


3.26E- 


08 


1.64E-08 


1.57E-08 


33 g 


1.71E-04 


1.11E-03 


8.41E-04 


1.31E-03 


66 Cu 


1.23E-10 


1.65E- 


10 


7.10E-11 


8.85E-11 


34g 


1.45E-04 


3.39E-05 


3.47E-05 


3.59E-04 


67 Cu 


1.21E-11 


2.06E- 


10 


6.61E-11 


7.34E-11 


35 g 


3.83E-10 


3.96E-10 


4.72E-10 


1.77E-09 


68 Cu 


9.76E-12 


1.65E- 


10 


5.32E-11 


6.91E-11 


36 g 


1.45E-09 


3.52E-10 


3.41E-10 


2.03E-10 


59 Zn 


1.07E-03 


3.81E- 


04 


6.15E-05 


6.95E-04 


37g 


1.87E-12 


2.00E-11 


1.48E-11 


3.27E-12 


60 Zn 


6.55E-01 


3.70E- 


01 


1.29E-01 


2.08E-01 


32 CI 


4.80E-04 


2.21E-05 


1.02E-05 


1.04E-05 


61 Zn 


1.35E-02 


2.74E- 


03 


1.51E-03 


1.78E-03 


33 CI 


9.07E-05 


6.11E-05 


9.32E-06 


7.32E-05 


62 Zn 


7.68E-01 


3.62E- 


02 


1.80E-02 


1.93E-02 


34 CI 


2.31E-04 


7.34E-05 


1.68E-05 


3.23E-05 


63 Zn 


2.37E-02 


1.06E- 


03 


6.10E-04 


5.23E-04 


35 CI 


3.10E-04 


1.45E-04 


5.12E-05 


2.93E-04 


64 Zn 


1.06E+00 


4.49E- 


02 


2.30E-02 


2.16E-02 


36 CI 


2.45E-07 


3.80E-08 


2.64E-08 


7.18E-08 


65 Zn 


2.37E-03 


1.00E- 


04 


5.13E-05 


4.82E-05 


37 C1 


1.46E-07 


1.60E-08 


1.35E-08 


3.73E-08 


66 Zn 


2.14E-03 


9.02E- 


05 


4.55E-05 


4.34E-05 


38 CI 


3.17E-12 


3.17E-11 


2.95E-11 


8.79E-12 


67 Zn 


6.02E-07 


2.57E- 


08 


1.30E-08 


1.24E-08 


34 Ar 


1.84E-03 


1.69E-04 


4.99E-05 


7.78E-05 


68 Zn 


1.09E-08 


6.34E- 


10 


3.16E-10 


3.11E-10 


35 Ar 


8.13E-04 


2.26E-04 


5.11E-05 


9.33E-05 


69 Zn 


9.27E-12 


2.04E- 


10 


6.15E-11 


7.40E-11 


36 Ar 


6.96E-02 


5.88E-01 


3.65E-01 


4.97E-01 


70 Zn 


1.00E-11 


2.42E- 


10 


6.58E-11 


8.77E-11 


37 Ar 


5.31E-05 


1.68E-04 


1.16E-04 


1.79E-04 


71 Zn 


6.19E-12 


1.08E- 


10 


4.16E-11 


5.50E-11 


38 Ar 


3.94E-04 


4.02E-05 


3.34E-05 


3.25E-04 


61 Ga 


8.83E-05 


3.80E- 


06 


1.89E-06 


1.84E-06 


39 Ar 


7.03E-08 


3.00E-09 


1.53E-09 


1.54E-09 


62 Ga 


4.83E-05 


5.15E- 


06 


1.67E-06 


4.16E-06 


40 Ar 


8.68E-08 


3.75E-09 


1.89E-09 


1.79E-09 


63 Ga 


1.37E-03 


1.39E- 


03 


4.86E-04 


9.56E-04 


41 Ar 


2.41E-12 


3.70E-11 


2.57E-11 


2.06E-11 


64 Ga 


9.88E-04 


4.41E- 


04 


3.13E-04 


2.06E-04 


42 Ar 


2.80E-12 


4.55E-11 


2.62E-11 


1.35E-11 


65 Ga 


7.21E-04 


3.52E- 


05 


2.43E-05 


1.63E-05 


43 Ar 


4.53E-12 


5.33E-11 


3.00E-11 


1.73E-11 


66 Ga 


8.58E-05 


3.66E- 


06 


2.18E-06 


1.77E-06 


36 K 


6.58E-04 


3.03E-05 


1.39E-05 


1.39E-05 


67 Ga 


5.36E-04 


2.26E- 


05 


1.15E-05 


1.09E-05 


37 K 


1.11E-04 


6.81E-05 


1.65E-05 


2.06E-05 


68 Ga 


1.44E-06 


6.12E- 


08 


3.08E-08 


2.95E-08 


38 K 


7.55E-05 


3.07E-05 


1.31E-05 


9.38E-06 


69 Ga 


1.02E-06 


4.32E- 


08 


2.18E-08 


2.08E-08 


39 K 


6.15E-04 


6.88E-05 


3.66E-05 


1.59E-04 


70 Ga 


6.56E-10 


2.11E- 


10 


9.27E-11 


1.24E-10 


40 K 


5.99E-07 


2.73E-08 


1.35E-08 


1.65E-08 


71 Ga 


1.99E-11 


1.74E- 


10 


7.51E-11 


1.13E-10 


41 K 


9.58E-07 


4.06E-08 


2.05E-08 


1.99E-08 


72 Ga 


9.93E-12 


1.38E- 


10 


5.53E-11 


7.30E-11 
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ivioaeis 


A 1 




R 1 

J3- 1 


R 




-D-O 






A 1 




J3- 1 


R 






42 K 


1.79E- 


11 


6.95E-11 


3.91E- 


11 


3.79E 


11 


73 Ca 


1.10E- 


11 


1.52E-10 


5.98E- 


11 


8.06E-11 


43 K 


7.21E- 


12 


8.11E-11 


4.19E- 


11 


3.41E- 


11 


63 Ge 


2.79E- 


05 


1.48E-06 


4.88E- 


07 


1.60E-06 


44 R 


4.45E- 


12 


7.90E-11 


3.84E- 


11 


3.77E- 


11 


64 Ge 


4.06E- 


02 


1.69E-02 


7.24E- 


03 


1.05E-02 


45 K 


4.16E- 


12 


1.00E-10 


3.78E- 


11 


4.05E- 


11 


65 Ge 


1.50E- 


04 


7.10E-05 


5.23E- 


05 


4.21E-05 


38 Ca 


1.42E- 


03 


9.17E-05 


3.61E- 


05 


3.80E- 


05 


66 Ge 


8.48E- 


03 


5.55E-04 


2.80E- 


04 


3.00E-04 


39 Ca 


3.00E- 


03 


2.06E-04 


7.89E- 


05 


8.39E- 


05 


67 Ge 


2.61E 


04 


1.18E-05 


6.26E- 


06 


5.58E-06 


40 Ca 


7.03E- 


02 


5.71E-01 


3.59E- 


01 


4.63E- 


01 


68 Ge 


4.68E- 


02 


1.98E-03 


9.97E- 


04 


9.52E-04 


41 Ca 


4.81E- 


05 


3.77E-05 


2.35E- 


05 


3.64E- 


05 


69 Ge 


2.27E- 


04 


9.59E-06 


4.84E- 


06 


4.62E-06 


42 Ca 


8.38E- 


04 


3.65E-05 


1.84E- 


05 


2.47E- 


05 


70 Ge 


1.04E- 


03 


4.38E-05 


2.21E- 


05 


2.11E-05 


43 Ca 


1.55E 


06 


6.57E-08 


3.37E- 


08 


3.32E- 


08 


71 Ge 


8.96E- 


07 


3.81E-08 


1.92E- 


08 


1.84E-08 


44 Ca 


1.14E- 


06 


4.81E-08 


2.42E- 


08 


2.33E- 


08 


72 Ge 


2.78E- 


08 


1.41E-09 


7.61E- 


10 


6.81E-10 


45 Ca 


6.55E- 


11 


1.45E-10 


5.04E- 


11 


7.25E- 


11 


73 Ge 


2.33E- 


11 


1.63E-10 


9.63E- 


11 


8.29E-11 


46 Ca 


5.62E- 


12 


9.28E-11 


5.45E- 


11 


4.61E- 


11 


74 Ge 


1.26E- 


11 


1.53E-10 


8.87E- 


11 


8.11E-11 


47 Ca 


4.46E- 


12 


6.63E-11 


3.61E- 


11 


4.57E- 


11 


75 Ge 


6.60E- 


12 


1.63E-10 


8.86E- 


11 


5.69E-11 


48 Ca 


4.01E- 


12 


8.17E-11 


4.34E- 


11 


3.33E- 


11 


65 As 


3.70E- 


09 


1.69E-10 


8.19E- 


11 


8.45E-11 


40 Sc 


8.98E- 


08 


3.80E-09 


1.91E- 


09 


1.83E- 


09 


66 As 


5.56E- 


07 


1.83E-07 


3.83E- 


08 


5.49E-08 


41 Sc 


1.51E- 


04 


6.64E-06 


3.22E- 


06 


3.16E- 


06 


67 As 


3.99E- 


05 


1.63E-05 


5.84E- 


06 


7.92E-06 


42 Sc 


2.13E- 


05 


1.05E-06 


4.67E- 


07 


4.56E- 


07 


68 As 


4.13E- 


06 


9.82E-07 


4.16E- 


07 


4.04E-07 


43 Sc 


1.43E- 


04 


7.64E-06 


3.29E- 


06 


3.55E- 


06 


69 As 


2.77E- 


06 


1.80E-07 


7.38E- 


08 


7.03E-08 


44 Sc 


1.56E- 


06 


6.62E-08 


3.34E- 


08 


3.25E- 


08 


7(1 As 


7.92E- 


07 


3.38E-08 


1.76E- 


08 


1.63E-08 


45 Sc 


8.43E- 


05 


3.56E-06 


1.87E- 


06 


1.72E- 


06 


71 As 


1.87E- 


05 


7.90E-07 


3.99E- 


07 


3.80E-07 


46 Sc 


3.40E- 


08 


1.53E-09 


7.62E- 


10 


7.26E- 


10 


72 As 


1.55E 


07 


6.74E-09 


3.38E- 


09 


3.25E-09 


47 Sc 


1.45E- 


09 


1.89E-10 


7.61E- 


11 


1.21E- 


10 


73 As 


1.09E- 


07 


4.92E-09 


2.52E- 


09 


2.37E-09 


48 Sc 


5.61E- 


12 


1.23E-10 


4.45E- 


11 


6.76E- 


11 


74 As 


1.50E- 


10 


1.38E-10 


8.36E- 


11 


9.52E-11 


49 Sc 


5.77E- 


12 


7.81E-11 


5.26E- 


11 


4.27E 


11 


75 As 


8.20E- 


11 


2.17E-10 


1.15E- 


10 


1.39E-10 


42 Ti 


5.85E- 


04 


2.46E-05 


1.22E- 


05 


1.18E- 


05 


76 As 


8.53E- 


12 


1.35E-10 


5.92E- 


11 


6.42E-11 


43 Ti 


1.23E- 


04 


6.28E-06 


2.55E- 


06 


2.96E- 


06 


67 Se 


2.21E- 


08 


9.79E-10 


2.79E- 


10 


3.95E-10 


44 Ti 


1.28E- 


02 


3.30E-03 


1.60E- 


03 


2.61E- 


03 


68 Se 


2.86E- 


04 


4.83E-05 


1.55E- 


05 


2.52E-05 


45 Ti 


7.11E- 


05 


2.69E-05 


1.52E- 


05 


1.20E- 


05 


69 Se 


1.26E- 


06 


3.68E-07 


1.39E- 


07 


2.36E-07 


46 T j 


1.16E- 


03 


6.18E-05 


3.25E- 


05 


3.07E- 


05 


70 Se 


1.24E- 


05 


7.52E-07 


3.38E- 


07 


4.00E-07 


47 T j 


4.09E 


06 


1.86E-07 


5.71E 


07 


9.57E- 


08 


71 Se 


1.00E 


06 


4.62E-08 


2.26E 


08 


2.16E-08 


48 Ti 


1.50E- 


05 


6.32E-07 


3.19E 


07 


3.05E 


07 


72 Se 


7.15E- 


04 


3.02E-05 


1.52E- 


05 


1.46E-05 


49 Ti 


1.43E- 


08 


7.22E-10 


3.61E 


10 


3.65E 


10 


73 Se 


8.71E- 


06 


3.68E-07 


1.86E- 


07 


1.77E-07 


50 T j 


3.54E 


11 


1.17E-10 


5.98E- 


11 


6.37E 


11 


74 Se 


1.48E- 


04 


6.27E-06 


3.16E- 


06 


3.02E-06 


51 Ti 


4.88E- 


12 


9.63E-11 


3.21E- 


11 


5.57E- 


11 


75 Sc 


9.18E 


07 


3.91E-08 


1.98E- 


08 


1.89E-08 


44 y 


3.86E 


04 


1.65E-05 


8.17E- 


06 


7.89E- 


06 


76 Se 


1.77E- 


06 


7.64E-08 


3.91E 


08 


3.69E-08 


45 y 


1.57E- 


04 


5.85E-05 


1.29E- 


05 


2.49E 


05 


77 Se 


2.21E- 


09 


2.86E-10 


1.28E- 


10 


1.08E-10 


46 y 


7.61E- 


05 


1.86E-05 


4.71E- 


06 


6.83E 


06 


78 Se 


1.88E 


10 


2.45E-10 


1.22E- 


10 


1.04E-10 


47y 


8.10E- 


04 


6.57E-05 


2.59E 


05 


3.11E- 


05 


69 Br 


4.76E- 


12 


2.07E-13 


1.03E- 


13 


1.01E-13 


48 y 


6.82E 


05 


2.94E-06 


1.54E 


06 


1.47E- 


06 


70 Br 


1.49E- 


09 


1.15E-10 


3.49E 


11 


4.13E-11 


49 y 


2.86E- 


04 


1.23E-05 


7.88E- 


06 


6.20E- 


06 


71 Br 


7.23E 


07 


4.28E-08 


1.86E- 


08 


1.93E-08 


50y 


9.02E 


08 


3.98E-09 


2.00E- 


09 


1.92E- 


09 


72 Br 


3.45E 


08 


3.46E-09 


1.46E- 


09 


1.89E-09 


Sly 


2.34E 


07 


1.00E-08 


5.03E 


09 


4.87E 


09 


73 Br 


5.42E 


08 


1.07E-08 


3.20E 


09 


5.88E-09 


52y 


2.79E- 


11 


8.86E-11 


5.14E- 


11 


6.41E- 


11 


74 Br 


9.04E 


09 


1.78E-09 


8.19E- 


10 


8.46E-10 


53y 


7.36E- 


12 


1.15E-10 


5.36E- 


11 


5.58E- 


11 


75 Br 


2.73E- 


07 


1.54E-08 


7.13E- 


09 


7.77E-09 


46 Cr 


7.13E 


04 


3.69E-05 


1.53E 


05 


1.68E 


05 


76 Br 


1.03E- 


07 


5.96E-09 


2.71E- 


09 


2.93E-09 
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Models 


A-l 


B-l 


B-2 


B-3 




A-l 


B-l 


B-2 


B-3 


47 Cr 

48 Cr 

49 Cr 


3.21E-04 
1.70E-02 
7.06E-04 


4.92E-05 
8.42E-03 
3.71E-04 


9.88E-06 
4.31E-03 
1.86E-04 


2.15E-05 
7.13E-03 
2.60E-04 


77 Br 

78 Br 

79 Br 


1.14E-06 
7.85E-09 
6.43E-09 


4.94E-08 
1.31E-09 
1.66E-09 


2.50E-08 
4.93E-10 
1.12E-09 


2.38E-08 
7.45E-10 
7.53E-10 



-35 - 



10.1 




8$ 1 ' 1 1 1 L 

J 4 5 6 7 B 



Logp c (gem 3 ) 

Fig. 1. — Evolutionary tracks of central temperature and density of the stars with 300M Q 
(thin dotted line), 5OOM (thick solid line), and 1000M© (thick dashed line). 
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Fig. 2. — Chemical composition just before the explosion (when the central density reaches 
10 10 g cm -3 ) of the 5OOM star (upper panel), and 1000M Q star (lower panel). The iron 
core occupies more than 20% of the total mass for both cases. 
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M r /M 

Fig. 3. — Density structure (top panel) and temperature structure (bottom panel) of 25M , 
500M Q , 1OOOM models. The horizontal axis is the mass fractions M r /M tota \. The vertical 
axis shows the density (top), and the temperature (bottom), respectively. The data of 25M 
is from Umeda & Nomoto (2003). 
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Fig. 4. — Models in which explosion occurs (filled circles) or not (filled squares), depending 
on two parameters 6j et and e for 1OOOM models. The other parameters are set at fi = lOe, 
Mbho = 1OOM , / = 0.01 (see Table 1, 2, and 5). 
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Fig. 5. — Snapshots of density structure showing how the jet is propagating at 30 second 
and 100 second after we started the calculation. The dial is normalized by the star's radius 
(~ 7.7 x 10 12 cm). Note that 0j et is set at 15° for both models, top left: 30s (A-l); top right: 
100s (A-l); bottom left: 30s (B-l); bottom right: 100s (B-l). 
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Fig. 6. — Initial radial positions of the matter that will be accreted into the central black 
hole for models A-l (left), A-2 (middle), and B-l (right). The blank region near the origin 
corresponds to the region of black hole initially formed. 




Log p |g cm" 3 ! Log p (g cm" 3 } Log p (g cm" 3 ) 



Fig. 7. — Maximum temperatures and densities of each mesh point for the directions of 
6 = 0°, 9 = 15°, 9 = 45° for models A-l (left), A-2 (middle), and B-l (right). 
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Fig. 8. — Distributions of elements: Fe-group elements for 8 = 0° (top left), a-elements for 
6 = 0° (top right), a-elements for 6 = 15° (bottom left), a-elements for 9 = 45° (bottom 
right), for Model A-l. Note that the mass range is set to 350 — 500M Q in the top right panel 
to see clearly the distributions of a-elements, while in others it is set to 100 — 600M© 



-42- 





Fig. 10. — Same as Figure 8, but for Model B-l. 
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Fig. 11. — Mass fractions of ejected a - elements and 56 Ni as a function of the direction 9. 
These figures are for models A-l, A-2, and B-l. These values are obtained by integrating 
over radial direction for each 9. 
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Fig. 12. — (a): [X/Fe] for Fe-group elements as a function of Y e for history A. (b), (c): Mass 
fractions of Fe-group elements, proton and 4 He as a function of Y e for history A. Both (b) 
and (c) are for history A. Some different elements are shown separately in (b) and (c) for 
clarity. 
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Fig. 13. — Abundance pattern of jet material averaged for 15 cases (5 Y e values times 3 
density-temperature histories) . 
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Fig. 14. — Total abundance patterns including jet contribution for higher resolution models. 
The open pentagons show the abundance ratios of gas of the central region in M82 (Ranalli 
et al. 2005). The bars show the range of abundance ratios observed in ICM (Baumgartner 
et al. 2004; Peterson et al. 2003). 




Fig. 15. — Same as Figure 14, but for lower resolutoin models. 
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Fig. 16. — Total abundance patterns including jet contribution for Case B. The bars show 
the observational ranges of EMP stars' abundances in Galactic halo (Cayrel et al. 2004). 



